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Interface exchange processes in LaAlOa/SrTiOa induced by oxygen vacancies 
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Understanding the role of defects in oxide heterostructures is crucial for future materials control 
and functionalization. We hence study the impact of oxygen vacancies (OVs) at variable concentra¬ 
tions on orbital- and spin exchange in the LaAlOs/SrTiOs interface by first principles many-body 
theory and real-space model-Hamiltonian techniques. Intricate interplay between Hubbard U and 
Hund’s coupling Jh for OV-induced correlated states is demonstrated. Orbital polarization towards 
an effective Cg state with predominant local antiferromagnetic alignment on Ti sites near OVs is 
contrasted with t 2 g{xy) states with ferromagnetic tendencies in the defect-free regions. Different 
magnetic phases are identified, giving rise to distinct net-moment behavior at low and high OV 
concentrations. This provides a theoretical basis for prospective tailored magnetism by defect ma¬ 
nipulation in oxide interfaces. 

PACS numbers: 73.20.-r,71.27.-|-a, 75.70.Cn 


I. INTRODUCTION 

There is strong evidence from recent experiments that 
oxygen vacancies play a crucial role in the physics of 
heterostructures between LaAlOs (LAO) and SrTiC^ 
(STO) [H-Q, as well as for surface and bulk joj 

features of pure STO. For instance, they may be rel¬ 
evant for ferromagnetic (FM) a nd superconducting or¬ 
der found in LAO/STO |10l - [l3| . OVs serve as electron 
dopants and can render an otherwise band-insulating en¬ 
vironment metallic [l3- The released charge from 0^“ 
fills the Ti(3d) shell, which has dP occupancy in stoichio¬ 
metric STO. This introduces effects of electron correla¬ 
tion in strontium titanate [T5l - [l^ . a compound adjacent 
to the Mott-insulating 3d^ RTiOs series (R: rare-earth 
ion) with a Ti(3d^) configuration. On similar grounds, 
the vacancy-induced doping enhances correlations in the 
LAO/STO interface. Puzzling interplay between itiner¬ 
ant and localized electrons in STO-based surfaces and 
interfaces is indeed suggested from scanning-tunneling 
spectroscopy [20l - [^ . magnetoresistance and anomalous 
Hall-effect measurements [2^, resonant x-ray scatter¬ 
ing [H, [2^ and photoemission 0, . 

A deeper comprehension of the defect influence on the 
interface phenomenology is motivated not just by basic 
research [27j |. Since the physical properties of oxide het¬ 
erostructures, ranging from insulating and/or conduct¬ 
ing to magnetic and/or superconducting, may be sub¬ 
tly tuned by the presence of impurities, promising en¬ 
gineering aspects emerge. Due to increasing control in 
detailed oxide-interface fabrication, new opportunities in 
high-response materials design are within reach [2^. In 
view of future spintronics devices, selective magnetic ac¬ 
tivation on the nano scale by versatile ways of defect 
creation [ 2 ^ may soon become available. 

Theoretical accounts of realistic LAO/STO het¬ 
erostructures are challenging because of the unique com¬ 
bination of complexity from the basic interacting quan¬ 
tum perspective and the structural bulk-to-interface set¬ 
ting. Calculations based on density functional theory 
(DFT) using hybrid functionals or employing static cor¬ 


relation effects from a Hartree-Fock-like treated Hubbard 
Hamiltonian (’+U’) can reveal some relevant aspects of 
the intriguing electronic structure (TB - OtL IsoI-I^. But 
there are two serious drawbacks to such extended Kohn- 
Sham schemes. Broken-symmetry states, i.e. long-range 
magnetic and/or charge orders, have to be stabilized of¬ 
ten right from the start to address correlation effects. 
Second, several many-body hallmarks such as param¬ 
agnetic local-moment behavior, low-energy quasiparticle 
(QP) formation, band narrowing, and interplay between 
QPs and Hubbard bands are not incorporated due to a 
lack of frequency dependence in the local electronic self¬ 
energy S. Invoking DFT-fdynamical mean-field theory 
(DMFT) overcomes these deficiencies, but, especially for 
defect environments treated by larger cells [ 13 , 1^, it 
remains numerically expensive. It is important to note 
that many-body physics beyond standard DFT-based ap¬ 
proaches is here not only a detail, but essential for illu¬ 
minating mechanisms of future technological use. As a 
further relevant aspect, first-principles supercell compu¬ 
tations to reveal defect influences are generally not per¬ 
fectly suited to the problem at hand. They are restricted 
by the choice of the (often too small) cell size and can 
introduce artifacts because of the introduced short-range 
defect ordering. 

In this theoretical work we want to focus on the pe¬ 
culiar problem of OVs at the LAO/STO interface over 
a larger concentration range. Though a deeper rele¬ 
vance of these defects for the formation of the original 
quasi-twodimensional electron liquid is still under de¬ 
bate, several first-principles calculations have shown that 
even the stoichiometric heterostructure is metallic from 
partly filled Ti(t 2 g)-dominated bands at the Fermi level. 
Although there are other theoretical suggestions [IB, isB - 
IBBl, OVs provide a natural way to explain the occurrence 
of interface ferromagnetism [TB, [H, HB, 113 ■ conven¬ 

tional DFT, an OV induces crystal-field lowered Cg-like 
impurity states on the neighboring Ti ions. Thinking in¬ 
tuitively, two limiting scenarios could apply depending 
on the concentration of vacancies in LAO/STO. Lin and 
Demkov [H, |4l| studied the dilute-defect limit with only 
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few oxygen defects, where Cg-like local moments on as¬ 
sumed Anderson/Kondo impurities may form. The lat¬ 
ter can couple ferromagnetically via Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interaction mediated by the itin¬ 
erant t 2 g electrons. On the other hand, in a dense- 
defect limit, the physics is closer to a minimal two-orbital 
(eg, t 2 g) Hubbard model near quarter filling [l^, Id^ . 
In a recent DFT-fDMFT work [l^ it was shown that 
in this limit, emerging FM order in the interface Ti02 
layer can indeed be explained by effective (Zener) double¬ 
exchange [H, 115 between an vacancy-induced e „ or bital 
and an in-plane t 2 g{xy) orbital. Michael! et al. also 
proposed Zener exchange between localized and itiner¬ 
ant states, but without referring to OVs as the source for 
localization. 

Although experimentally the influence of OVs in STO- 
based materials is documented by monitoring pWsical 
response with varying oxygen partial pressure [J, lly , 
a good quantitative understanding of the OV concentra¬ 
tion as well as definite information on the location with 
respect to the interface ii is still lacking. To cope 
with the uncertainties in the number of OVs, we here 
perform investigations in a broad concentration range. 
In order to achieve this task, the correlated electronic 
structure is treated in a real-space framework allowing 
for in principle arbitrary vacancy configurations in num¬ 
ber and arrangement. Depending on the concentration of 
OVs, we encounter different orbital- and spin-exchange 
regimes, that shed light on the emerging FM order at 
the LAO/STO interface. Local and non-local processes 
subject to a subtle interplay between the Hubbard U and 
Hund’s Jh govern the OV-induced magnetism. Through¬ 
out the paper all energies are given in electron volts. 

The paper is organized as follows. To set the stage, we 
touch base with previous DFT-t-DMFT calculations (l9j | 
and start in Sec. HIl with a brief view on the impact of the 
Hund’s coupling Jh on the magnetic order in a n-type 
LAO/STO interface in the limit of high OV concentra¬ 
tion. In Sec. IIII Al the correlated real-space modeling 
to describe different defect concentrations is introduced. 
Results in the dilute limit of a single oxygen vacancy 
as well as of two OVs in the interface are presented in 
Sec. IIII HI Section fill Cl deals with the evolution of elec¬ 
tron correlation and magnetism upon increasing the num¬ 
ber of OVs from the dilute- to a dense-defect limit. The 
work closes with a discussion and summary in Sec. IIVI 


II. INFLUENCE OF Jh IN THE DENSE-DEFECT 
LIMIT OF OXYGEN VACANCIES IN LAO/STO 

A charge self-consistent DFT-I-DMFT study in a 
dense-defect limit of 25% OVs in the Ti02 interface layer 
was performed in Ref. [l^. We define that limit by OVs 
exclusively located in the Ti 02 interface layer, with each 
Ti ion having one OV in bonding distance (see Fig. [TJ. 
Here and throughout this paper Ti neighborhoods with 
more than one nearby OV are not considered. 



FIG. 1. (color online) (a) Ti02 interface layer of n-type 
LAO/STO with 25% oxygen vacancies utilized in supercell 
DFT-fDMFT calculations [I^. Ions Til and Ti2 form the 
basis in the V2xy/2 primitive cell, (b) minimal relevant Ti 
orbitals with |eg} ~ 0.55|z^) ± 0.84|a;^—i/^} [l9l| . 


For the DFT part a mixed-basis pseudopotential 
framework is used and the DMFT impurity problems are 
solved by continuous-time quantum Monte Carlo (CT- 
QMC) A minimal correlated subspace was de¬ 

rived to consist of a two-orbital [e~g, t 2 g{xy)] manifold 
located at the interface Ti ions. Remaining t 2 g orbital 
degrees of freedom are included in more distant lay¬ 
ers from the interface. The local Coulomb interactions 
in the interacting Hamiltonian with Slater-Kanamori 
parametrization, i.e., the Hubbard U and Hund’s ex¬ 
change Jh were set to U=2.5 and Jh=0.5, in line with 
other works M- A double-exchange-like (DE) mecha¬ 
nism is effective in stabilizing FM order within the in¬ 
terface. Thereby the spin polarization is triggered by 
inter-orbital scattering between the nearly-localized e~g 
state and the less-localized xy state. Because of this 
exchange mechanism, the fewer and more-itinerant xy 
electrons exhibit stronger spin polarization. The self- 
consistently adapting Ti occupation results in a quarter 
filling (nTii,Ti 2 "^l) of the two-orbital correlated subspace 
within the interface layer. In total, each OV releases 2 
electrons, which add to the 0.5 electrons per Ti ion due to 
the polar-catastrophe avoidance within LAO/STO het¬ 
erostructures. In the supercell treatment of the dense- 
defect limit, which notably includes several Ti 02 layers, 
these two contributions adjust such that I electron per Ti 
settles in the interface Ti02 layer. Hence the latter layer 
formally consists of Ti^+(3d^) ions in the given defect 
limit. 

Since the size of the Hund’s coupling is key to a 
double-exchange-like mechanism for the onset of ferro¬ 
magnetism, we here present results from varying Jh, 
while keeping U=2.5 fixed. For more details on the cal- 
culational setup see Ref. [H. Usually the value of Jh 
is much less modified by screening processes than the 
Hubbard U and therefore remains close to the magni¬ 
tude in the free atom. Hence although for LAO/STO 
a Hund’s coupling Jh=0.5— 0.6 is expected, we take the 
liberty of changing this value in order to assess it rele¬ 
vance for the underlying physics. Figure[2Ka) displays the 
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the FM order for Jh=0.7 is dominantly carried by the 
xy orbital. These findings underline the importance of 
DE processes in the formation of ferromagnetism at the 
LAO/STO interface in the dense-defect limit. 


III. REAL-SPACE APPROACH TO OXYGEN 
VACANCIES IN LAO/STO 


The scope is now broadened by addressing lower OV 
concentrations in the Ti02 interface layer. Performing 
large-scale incommensurate studies from first principles 
of the interplay between defects, disorder and correla¬ 
tions for multi-orbital lattice is nowadays still too expen¬ 
sive. We thus need a minimal setting that is geared to 
carry the key physics of the OV-doped LAO/STO inter¬ 
face. Therefore a model-Hamiltonian approach is used, 
inspired by the results of the DFT-l-DMFT investigation 
in the previous section. 


FIG. 2. (color online) DFT-DMFT results depending on Jh 
for the dense-defect LAO/STO interface with magnetic order 
(r=180K). (a) Orbital occupations (left) and Ti magnetic 
moment (right), (b) Total (left) and local (right) spectral 
functions. 


orbital- and spin-dependent occupations with Jp in the 
magnetically ordered phase of the dense-defect limit. For 
Jh^O. 4 ferromagnetism with local Ti moment m~0.2- 
0.4 /xb occurs. But when Jp becomes smaller and even¬ 
tually tends to zero, antiferromagnetic (AFM) order be¬ 
tween the nearest-neighbor (NN) titanium ions sets in. In 
the latter regime the orbital polarization in favor of the 
more localized e'g level strongly increases towards nearly 
full polarization at Jp=0. 

Figure [2Kb) shows the fc-integrated spectral function 
Acr(w)=^jj. Acr(k, w) in the limiting cases Jr= 0, 0.7 for 
the whole supercell as well as for the correlated subspace 
of the Ti[e'g, t2g{xy)] states. Concerning the correlation 
strength the Hund’s coupling has a known model impact 
within a two-orbital system near quarter filling (50j |. In 
the case of vanishing Jp stronger correlations occur, giv¬ 
ing rise to a prominent lower Hubbard peak at 1.3. 
This incoherent excitation is exclusively associated with 
the vacancy-induced e~g state and resembles a similar fea¬ 
ture in photoemission data 0, H, [1^ . Locally, hopping 
is nearly blocked on the interface Ti ions and residual 
metallicity is mainly provided by sites far from the inter¬ 
face. For rather large Jp the interface is well conductive, 
with a prominent QP peak at the Fermi level but an ab¬ 
sent lower Hubbard peak. Coexistence of general metal¬ 
licity with the latter incoherent excitation holds for rea¬ 
sonable intermediate values of Jp M- The AFM order 
for Jp=0 is completely carried by the e"g orbital, while 


A. Model and methodology 


A two-orbital Hubbard Hamiltonian based on the 
vacancy-induced effective e~g state and the in-plane 
t 2 g{xy) state is employed on a 10x10 square lattice with 
fVTi=100 titanium ions, mimicking the interface Ti 02 
layer (see Fig. IH). Only the Ti sublattice is treated ex¬ 
plicitly and the oxygen degrees of freedom are integrated 
out within the chosen Hamiltonian form. Since the cre¬ 
ation of OVs amounts to electron doping, explicit involv- 
ment of remaining oxygen orbitals can be neglected to 
first approximation. Periodic boundary conditions are 
applied. Only intra-orbital NN hoppings are retained in 
the model. 


Lets focus first on the lattice scenario in the dense- 
defect limit, where each Ti site is affected by a nearby 
OV, to build up the model characteristics. In line 
with Ref [l^, we choose te =txy=0.2 for the NN hop¬ 
pings. In contrast to a different modeling by Pavlenko 
et al. [1^, our hoppingamplitudes from the projected- 
local-orbitals method ^ for higher OV concentrations 
are not strongly orbital dependent. From a noninteract¬ 
ing point of view, the crystal-field splitting A between the 
xy level and the vacancy-induced low-energy e~g is the key 
model parameter. Note that A is different from the usual 
octahedral crystal-field splitting that is already vital in 
the stoichiometric compound. The latter energy splitting 
does not occur in the present defect model. Again from 
Ref we set A=0.3. 


Local electron-electron interactions via the Hubbard U 
and Hund’s coupling Jp are applied at all Ti sites. The 
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full Hamiltonian then reads 


H — t ni^xy,a) + 

(•i,j)acr icr 

^ ^ “t“ 


E {^' '^iota'^ioL' (7 + u” (1) 

i,a^a' ^cr 

^iacr^iaa^ia'a^ia' J* ’ 


where c(c^) are creation (annihilation) operators, n=c^c, 
i, j are site indices, a, a'=l3,xy and cr=t, i marks the spin 
projection, using U'=U — 2Jh, U"=U — 3Jh- For the 
same Hubbard H, the strength of electronic correlations 
is usually weaker within slave-boson theory than within 
CT-QMC. If not otherwise stated the Hubbard U is thus 
set to U=i in all calculations. The Hund’s coupling is 
set to Jh=0.55, again with variations to smaller/larger 
values to trace its relevance. For the dense-defect limit, 
Ai=A and /3=eg, i.e. the vacancy-induced Cg crystal- 
field state is active on each Ti site. 

In the latter limit every interface Ti ion has one neigh¬ 
boring OV and the given Hamiltonian is coherently appli¬ 
cable on each Ti site. However, when the defect number 
is reduced, titanium sites without a nearby OV appear 
and at those Ti sites there are no low-energy eg states. 
In the stoichiometric case the eg orbitals are strongly 
bound to 0(2p) and do not contribute either to states 
at the Fermi level or to any possible local-moment for¬ 
mation. In order to keep the modeling simple, we make 
the following approximations when treating general de¬ 
fect cases in real space: (i) the Hamiltonian of form o 
is used throughout the lattice, (ii) the parametrization of 
Ai is 


. _ J A , /3 = Cg : if an OV nearby , . 

* 0 , /3 = xz/yz : if no OV nearby ' ' 

and (iii) multiple OVs around a Ti site are prohibited. 
The interpretation of (ii) is as follows: without a nearby 
OV, the Ti local low-energy electrons are mainly of t 2 g 
kind and thus the former e~g degree of freedom takes over 
the role of an additional effective t 2 g orbital. This can 
be justified by a notable hybridization between e"g and 
xz,yz in the dense-defect case [l^. We neither change 
hoppings for Ti sites with or without nearby OVs nor 
employ a concentration-dependent hopping modification. 
Such a more detailed parametrization is hard to fix and 
the aim here is to work in a canonical two-orbital setting. 

The final modeling step examines the concentration- 
dependent electron filling of a lattice with vacancy con¬ 
centration c=N^ac/No, where Wac is the number of OVs 
and Ao=2VTi denotes the number of oxygen sites. Our 
electron count considers only the single Ti 02 interface 
layer, explicit charge fluctuations to more distant layers 
across the interface are neglected. In the dense-defect 
limit of the supercell treatment, DFT-I-DMFT yields a 



FIG. 3. (color online) Real-space two-orbital modeling on a 
10 X10 Ti02 square lattice for the n-type LAO/STO interface. 


filling of one electron at each Ti site in the interface 
Ti02 layer. Supercell DFT calculations for the defect- 
free interface result in an itinerant two-dimensional elec- 
tro n sy stem with a count of 0.5 electrons per interface 
Ti [^, i.e. 50 electrons for the 100 Ti sites, in line with 
the polar-catastrophe avoidance. Putting these num¬ 
bers together, we choose the linear-interpolation scheme 
»^tot=.^Ti/2-l-iVvac for the total lattice electron count ritot 
at intermediate defect levels. 

The interacting multi-orbital problem is solved by a 
real-space formulation of the rotational-invariant slave- 
boson (RISE) mean-held method [13, It corre¬ 

sponds to single-site DMFT close to zero temperature 
with a simpler impurity solver than the CT-QMC, al¬ 
lowing for local self-energies S with a linear frequency 
dependence and static terms. Renormalized QPs as well 
as local multiplets can be monitored in the interacting 
regime. Explicit intersite self-energy terms are neglected, 
but due to the coupling of all sites in the RISE self- 
consistency cycle effects of incoherency due to the distri¬ 
bution of defects are included. Our real-space approach 
is reminiscent of a single-orbital variant put into prac¬ 
tice by Andrade et al. [s^. However instead of simpli- 
Hed Kotliar-Ruckenstein slave bosons [13 we here use 
the full rotational-invariant extension and elaborate on 
a multiorbital framework. Eecause of the model/method 
complexity no disorder averages are performed in this 
work. The calculations utilize 60 slave bosons and 17 la- 
grangian multipliers per site, 7700 variational parameters 
in total, with dimension 400x400 for the kinetic part of 
the Hamiltonian. 

In the following, the ordered magnetic moment m, the 
orbital moment r, the paramagnetic local spin moment 
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mpM and the orbital polarization are defined as 
m = '^ma = - rial) > - nxy,a) 



where S denotes the local spin operator and 0={0). 
Lattice-averaged values Qiat of these quantities Q are 
computed by Q\e,t=Q 


B. Dilute-defect limit 

Since our parametrization is adjusted to the 
DFT-I-DMFT results within the dense-defect limit, the 
model performance in the contrary limit of one or two 
OVs is of primary interest. For the case of two vacan¬ 
cies, a long-distance accommodation is chosen. Further¬ 
more, they are placed on differently directed bars on the 
effective square lattice to minimize coherency effects. 

Figure |4] summarizes the obtained real-space results 
for key quantities in the PM phase as well as with mag¬ 
netic order. Each small square resembles one Ti site and 
the oxygen sites are located in the middle of the enclosing 
bars. As expected, the orbital polarization towards the e~g 
level for the Ti sites near vacancies is substantial. Away 
from the defect, especially in the two-OV case, there are, 
however, also oscillations in the orbital moment. The 
local PM spin moment is largest close to the OVs, un¬ 
derstood from the well-localized correlated eg electrons, 
but attains a minimum just at second-nearest Ti sites. 

When initializing the RISE calculations with FM or¬ 
der, at saddle-point the solutions indeed reveal a small 
net FM lattice moment. But throughout the lattice the 
site-resolved ordered magnetic moment m alternates. In¬ 
terestingly, while m for the single-OV case is still rather 
small at the defect, already the two-OV case displays a 
strong increase in the absolute value of the near-defect 
Ti magnetic moment. Furthermore, whereas for the sin¬ 
gle vacancy the spins on both NN Ti sites align in the 
same direction, an AFM alignment emerges with two va¬ 
cancies. This is surprising since the nearby Ti sites in 
both scenarios show similar strong filling of the localized 
eg state. Due to the strongly correlated regime, a signif¬ 
icant kinetic exchange /U of AFM type is expected. 
Therefore the distance between OVs has to matter for 
the near-defect spin coupling. Structures with defects at 
closer range favor the local AFM alignment. The non¬ 
local exchange with the alternating m throughout the 
lattice may also play a role in the local spin alignment. 
Thus there appears to be an intricate coupling between 
short- and longer-range exchange on the lattice. 

Note that in general the overall system is rather sus¬ 
ceptible to net magnetic order even with only a small 
number of defects and small resulting moments. The 
electron count per site reads h=0.51(0.52) for one(two) 


OV(s), i.e. in principle the two-orbital Hubbard model¬ 
ing is close to one-eighth filling. Thus though the stoi¬ 
chiometric interface in DFT-l-DMFT does not show mag¬ 
netic order [I^, it apparently only needs a low concen¬ 
tration of OVs to create the first magnetic instability in 
the correlated regime. 

In order to investigate the issue of magnetic exchange 
with distance in more detail, the averaged ordered mag¬ 
netic moment with spacing r from the OV in the single¬ 
defect case is plotted in Fig.O The NN Ti site is located 
in distance rNN=0.5 in units of the lattice constant a. 
Though it is expected that results will dependent on the 
chosen lattice size, there are significant variations within 
the shorter- and longer-range regions. For d^2.5 the 
moments predominantly change their sign and the spin 
coupling switches from FM- to AFM-like. Let’s assume 
here a Fermi-wave-vector kp modulated RKKY exchange 
with J{d)^cos{2kp d) in the low-density limit for the sea 
of conduction electrons [s^. This would correspond to 
kp^n/lO in reciprocal units, more or less commensurable 
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FIG. 4. (color online) Site-resolved quantities of interest for 
one OV (left) and two OVs (right) in the Ti02 layer, (a) PM 
orbital moment r, (b) PM local spin moment mpM, and (c) 
ordered magnetic moment m. 
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FIG. 5. (color online) Distance-dependent magnetic moment 
m{r) with the oxygen vacancy at the origin in the magnet¬ 
ically ordered single-defect case for different choices of the 
Hund’s coupling Jh. 


with the noted one-eighth filling. But one has to keep in 
mind that the local ordered moment in our single-OV 
case is too small for a serious application of the standard 
RKKY picture. Still, for the record, there are indications 
of RKKY-like exchange taking place in the dilute-defect 
limit. Concerning the influence of Hund’s coupling, not 
surprisingly, a larger Jh increases the spin moments since 
it tends to locally align the contributions from the orbital 
degrees of freedom. 


C. Electron correlation and magnetism from the 
dilute- to the dense-defect limit 

We now investigate the regime between the dilute- 
defect limit (c=0.005, 0.01) and the dense-defect limit 
(c=0.25). Additional configurations with randomly dis¬ 
tributed OVs for intermediate concentrations are con¬ 
structed, prohibiting local Ti neighborhoods with more 
than one vacancy in NN distance, respectively. In the fol¬ 
lowing, site-resolved and site-averaged data is presented. 
It proves instructive not only to perform full lattice aver¬ 
ages, but also to differentiate between the two groups of 
Ti sites, i.e., those with and those without nearby OVs. 

Figure [6] provides a measure of the correlation 
strength by displaying the paramagnetic QP weight 
~ Because of the low electron count 

at low OV concentrations the lattice QP weight starts 
off with values close to the noninteracting limit Ziat=l. 
With more vacancies and increased electron doping, gen¬ 
eral electronic correlations become stronger down to 
^iat~0.4 in the dense-defect case. For the considered 
electron-doping regimes the Hund’s coupling Jh on aver¬ 
age weakens correlations, in line with DFT-t-DMFT re¬ 
sults from Sec. im For any doping, orbital- and/or site- 
selective Mott transitions remain absent. Not surpris¬ 


ingly, electrons in the Cg orbitals are nonetheless stronger 
correlated since they are more localized due to the low¬ 
ered crystal field. While for the xy state an obvious site 
discrimination occurs only at higher doping, the Cg elec¬ 
trons near OVs are already much heavier at the lowest 
doping. Interestingly the xy electrons eventually become 
stronger correlated away from near-OV titanium. Thus 
the vacancy influence on correlations is twofold: it locally 
strengthens the Cg correlation and nonlocally fosters the 
xy correlations. 

As mentioned, the site-dependent correlation strength 
is of course related to the local orbital filling, shown in 
Fig. [3 from low to higher OV-induced electron doping. 
For comparison, orbital occupations in the noninteract¬ 
ing case (C/=Jh=0) are additionally depicted. For the Ti 
sites with nearby OVs a strong orbital polarization to- 



c 



FIG. 6. (color online) Orbital-resolved QP weight Z vs con¬ 
centration of oxygen vacancies (OVs) in the PM phase. Total 
electron count increases from ntot=51 (1 OV, c=0.005) to 
ntot = 100 (50 OVs, c=0.25). (a) Average Za for states on Ti 
sites with and without nearby OVs for Jh=0.55 (left) and 
Jh=0.7 (right), (b) Site-resolved Z for selected dopings and 
Jh=0.55 (top) as well as lattice Ziat for different choices of 
Jh (bottom). 
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FIG. 7. (color online) Orbital-resolved occupation in the PM 
phase, (a) Average occupations of Ti sites with and with¬ 
out nearby OVs for the noninteracting (left) and interacting 
(right) case, (b) Concentration-dependent orbital moment r 
in real space (top) and lattice-averaged (bottom). 


FIG. 8. (color online) Orbital-resolved magnetic moment m. 
(a) Average magnetic moment of Ti sites with and with¬ 
out nearby OVs for Jh=0.55 (left) and Jh=0.7 (right), (b) 
Goncentration-dependent magnetic moment m in real space 
(top) and lattice-averaged (bottom). 


wards Cg already at low dopings is obvious, even without 
local Coulomb interactions. Rather independent of the 
number of vacancies, a polarization C'^2.8 holds for these 
sites in the noninteracting problem. With interactions 
this orbital polarization is substantially increased for all 
doping levels because of the crystal-field renormalization 
via the electronic self-energy. It still gradually decreases 
from ^^6 in the dilute limit to ^^^4 in the dense limit. 

Due to the absence of local crystal fields a subtle com¬ 
petition between both orbital degrees of freedom occurs 
at the remaining Ti sites. However, remember that in 
this region the ’e~g’ orbital inherits the role of an addi¬ 
tional t 2 g orbital in our modeling. The noninteracting 
case does not reveal any finite orbital moment for any 
OV concentration. By including Coulomb interactions, 
low doping from the dilute limit slightly disfavors the 
xy orbital. But interestingly, at Cp~0.08 the xy orbital 
takes over the lead in occupation. Thus also here a non¬ 
local impact of the OV-induced correlations shows up, 
breaking the orbital degeneracy at Ti sites away from 
the defects. Since in the dense-defect limit every Ti ion 
has one nearby OV, the class of defect-free Ti sites dis¬ 
appears and its nominal filling then, of course, vanishes. 


Intuitively, the doping-dependent change of orbital-filling 
hierarchy for this Ti class can be understood as follows. 
At very low dopings, the interacting system tries to put 
more electrons into the e~g level, since there they can oc¬ 
casionally enjoy the lower crystal field close to an OV. 
Yet at some concentration with increased electron filling, 
it is more beneficial to put electrons which like to visit 
defect-free regions in the overall less occupied xy levels 
to minimize the Coulomb interaction and gain kinetic 
energy. 

The real-space variation of the site-dependent orbital 
moment r underlines these findings [see Fig. [71[b)]. Be¬ 
yond the critical doping level Cp there is a qualitative 
change in the polarization of the ’interstitial’ region to¬ 
wards xy. A shoulder in the lattice-averaged orbital 
moment riat is located around Cp. Of course, every in¬ 
crease in OVs, with strong orbital polarization towards Cg 
nearby, renders riat further monotonically growing with 
c. Give or take, the averaged influence of Jr on the or¬ 
bital moment is as expected from multiorbital Hubbard 
models, i.e. it works against the crystal field and tries to 
wash out orbital polarization [ssj) . 

The orbital- and site-resolved magnetic moment ex- 
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FIG. 9. (color online) Low-energy PM spectral informa¬ 
tion for dopings c=(0.045, 0.09, 0.175) (left to right) using 
Jh=0.55. Orbital-resolved site-contribution to the 

renormalized QP state in an energy 

window [—0.05, 0.05] around the Fermi level. 


hibits an even more intricate structure with respect to 
the vacancy concentration (cf. Fig. |S]). Starting from 
the dilute-defect limit the Ti sites with nearby OVs can 
be divided into two subclasses. One spin thereof compen¬ 
sates their net moment by AFM alignment, the other one 
displays an FM alignment with smaller moments. There¬ 
from a still sizable net FM moment near c^O.05 results, 
whereby the RKKY-like exchange noted in Sec. lIIIBl can 
be held responsible. Intriguingly, close to Cp this net FM 
moment vanishes, accompanied by the disappearance of 
the Ti-site subclass with FM alignment near OVs. In ad¬ 
dition, the nearly exclusive AFM alignment at OVs above 
Cp comes along with zero spin polarization on the re¬ 
maining Ti sites in the concentration range 0.08<c<0.13. 
Thus the RKKY-like driven FM phase is followed by a 
phase region of separated AFM pairs with zero net mo¬ 
ment. Though, as discussed before, the xy orbital polar¬ 
ization in the ’interstitial’ is already active, in this con¬ 
centration range there is no novel exchange mechanism 
yet strong enough to spin polarize the Ti sites aways from 
the defects. This changes for cde^O. 13, when eventu¬ 
ally the latter sites build up a considerable a:?/-dominated 
magnetic moment. Effective non-local double-exchange 
between the more-localized e'g electrons near OVs and 
the more-itinerant xy electrons away from OVs yields 
net ferromagnetism. Also, part of the AFM alignment of 
the near-defect Ti sites breaks up and these sites join in 
the contribution to the FM order. Interestingly, there is 
a minimum in the averaged eg magnetic moment on the 
way to the dense-defect limit. This may be explained by 


the competition between FM double-exchange and AFM 
kinetic exchange on Ti near OVs. The increasing elec¬ 
tron filling re-strengthens the kinetic exchange close to 
the defect in the now FM-polarized environment, i.e. a 
re-formation of AFM pairs occurs. A larger Hund’s cou¬ 
pling shifts that minimum to higher electron dopings [see 
Fig. EKa)], i.e. the local Cg occupation has to be closer 
to the kinetic-exchange favored half-filling [cf. Fig. |71[a)] 
to overcome the Jn-supported DE processes. We also 
checked the generic influence of a smaller U=2.5 and 
encountered an overall somewhat reduced magnetic mo¬ 
ment and a weakening of the xy magnetism based on the 
nonlocal double-exchange. 

For the FM phases, our revealed lattice magnetic mo¬ 
ment 77iiat^0.1-0.2/rB is in very good agreement with ex¬ 
perimental findings 1^ . So for moderate Jr the nonlocal 
polarization effect of OVs fosters a sizable xy-dominated 
magnetic moment on defect-free Ti sites above a concen¬ 
tration cde~0.I3. This nonlocal double-exchange pro¬ 
cess extends the local DE mechanism from coherent sys¬ 
tems, here active in the dense-defect limit. Note that 
the obtained magnetic moment m~0.2/rB for Jh= 0.55 in 
the latter limit is in excellent agreement with the for¬ 
mer DFT-fDMFT results (cf. Fig. [2]), highlighting the 
consistency of the model. However, there is a difference 
in the orbital contributions, since in the real-space RISE 
model the e~g level is more strongly spin-polarized than 
xy. This may be explained by the fact that fluctuations 
and their correlations, relevant for assessing the DE pro¬ 
cesses in detail, are underestimated in simplified RISE 
compared to DMFT with a CT-QMC solver. It could 
also be that scattering in additional Ti 02 layers, which 
is not included in the real-space modeling, supports the 
xy spin polarization. 


Finally, we address spectral features at low energy 
since, e.g., the question arises about the different site 
and orbital contributions to the resulting metallicity. For 
selective dopings. Fig. [H] shows the total QP density of 
states (DOS) as well as the site- and orbital-resolved QP 
spectral weight within a small energy window around 
the Fermi level, both in the PM regime. The real- 
space/orbital resolution is naturally derived from ana¬ 
lyzing the low-energy lattice eigenvectors of the renor¬ 
malized kinetic Hamiltonian. 


At the concentration just above Cp the total DOS dis¬ 
plays a pseudogap-like feature at the Fermi level. From a 
Stoner argumentation, this reduced spectral weight at £f 
is in line with a vanishing of ferromagnetism in this con¬ 
centration regime. In the case of higher OV numbers, i.e. 
higher electron dopings, the low-energy density of states 
rises again within the DE-FM region. As generally ex¬ 
pected the xy weight is mainly located in the defect-free 
regions, and the e~g weight stems dominantly from defect- 
near regions. Close to the dilute-defect limit, the overall 
xy low-energy weight is lower, but it becomes dominant 
just above Cp. There the e'g electrons appear most lo¬ 
calized, giving eventually rise to the pair-AFM phase. 
For cde^O.13 eventually the e'g low-energy contribution 
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again overcomes the xy one, marking the intriguing scat¬ 
tering regime of the effective double-exchange region. 


IV. SUMMARY AND DISCUSSION 

This work examined the key effects of the electronic 
structure reconstruction in the LAO/STO interface due 
to the presence of OVs. Different orbital and spin ex¬ 
change processes are identified for varying OV concen¬ 
trations. From the revealed and expected magnitudes 
of hoppings, crystal fields and Coulomb interactions a 
straightforward picture of fully localized (Kondo-like) 
electrons near OVs is not evident. Charge self-consistent 
DFT-I-DMFT for LAO/STO supercells in a dense-defect 
limit yields a dichotomy between OV-induced heavier e/ 
states as well as xy states with a higher QP weight. But 
itinerancy remains a common feature when including a 
finite Hund’s coupling Jr. The latter not only triggers 
the correlation strength but fosters ferromagnetism above 
a threshold doping through effective nonlocal and local 
double-exchange processes. Additionally, the strongly 
correlated dense-defect limit is spectrally marked by a 
lower Hubbard peak of e~g kind. For generic OV num¬ 
bers in a 10 X10 Ti 02 model interface our correlated real- 
space RISE approach elucidates further mechanisms due 
to the system separation into Ti sites near and away from 
OVs. Coulomb interactions are relevant to trigger intri¬ 
cate (nonlocal) orbital polarization processes, especially 
at the defect-free Ti sites. 

Concerning magnetic order, the schematic phase dia¬ 
gram in Fig. [To] summarizes the main findings. Already 
in the dilute limit of very few vacancies, oscillations in 
the sign of the magnetic moments with distance from the 
OV can be detected. Near OVs two subclasses of Ti sites 
appear, one favors local AFM and the other local FM 
alignment. An effective RKKY-like exchange mechanism 
weakly spin polarizes the system, giving rise to a finite 
FM net moment. We note that the present RKKY or¬ 
dering is not conventional in the sense that the involved 
local moments are comparatively small. Due to the del¬ 
icate exchange mechanism, the expected Curie tempera¬ 
ture Tc associated with this phase is rather low. Above 
an interface vacancy concentration Cp^O.08, the pairs 
of AFM-aligned Ti sites dominate the lattice and the 
spin polarization in the regions without OVs disappears 
until at cde~0.13 the double-exchange becomes strong 
enough to polarize the ’interstitial’, now with dominant 
xy character. In addition, the DE mechanism is effec¬ 
tive in switching local AFM pairs to FM alignment. A 
more robust ferromagnetic order sets in, with a suppos¬ 
edly much larger and continues to be stable up to the 
dense-defect limit. The value of Jr influences the compe¬ 
tition between eg filling-controlled re-strengthened AFM- 
like kinetic exchange and the DE processes near OVs in 
this novel DE-FM phase. We did not delve into the pos¬ 
sible phase transitions among the three phases. Because 
of the overall defect system without a coherent local or- 



FIG. 10. (color online) Rough sketch of a finite-temperature 
(T) magnetic phase diagram for OV concentrations c based 
on the real-space RISE modeling of the LAO/STO interface. 
Dashed phase boundaries of the pair-AFM phase mark the 
possible occurrence of coexistence regions with the nearby 
FM phases. 

der parameter and straightforward symmetry distinction, 
first-order (-like) transitions with coexistence regions are 
expected (cf. Fig. ITO)) . 

Our obtained behavior with increasing OV concentra¬ 
tion shares several features with experimental results. It 
is generally in accordance with the found key dependence 
of magnetism on electron doping in STO-based materi¬ 
als. An interplay of AFM and FM tendencies has been 
recently identified by Bi et al. [1^. The different ex¬ 
perimental results concerning the range of stability for 
LAO/STO ferromagnetism may be related to substantial 
differences in the number of vacancies in the respective 
samples. Whereas nearly stoichiometric interfaces are 
susceptible to the low-Tc RKKY-like FM phase [5^. 1^. 
OV-rich samples can stabilize the DE-FM phase with the 
surprisingly high near room temperature fl^ . . Of 
course, such different ferromagnetism may also emerge 
in very inhomogeneous samples. A theoretical resolu¬ 
tion of phase separation on a larger lattice scale based 
on the present modeling is, however, numerically hard 
to achieve. Concerning the crucial concentration depen¬ 
dence, unique behavior connected to critical electron den¬ 
sities has been revealed in magnetotransport measure¬ 
ments [ 2 ^. In that respect it would be very interesting 
to trace in detail the ferromagnetism in applied magnetic 
field B within a group of samples with different OV con¬ 
centrations, or to perform in-situ monitoring with oxy¬ 
gen pressure. For instance, the pair-AFM phase could be 
transformed to FM order by a larger field B. 

We have shown that itinerancy, polarized orbital de¬ 
grees of freedom and magnetic order naturally go to¬ 
gether in LAO/STO interfaces with OVs. The realistic 
many-body physics remains challenging and needs fur¬ 
ther work. From theory, the treatment of the full Ti(3d) 
shell for the demanding supercell- and real-space compu- 
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tations would allow for a more detailed orbital resolution. 
Including the impact of spin-orbit coupling and the com¬ 
petition of the revealed phases with superconductivity is 
a natural further modeling step. In addition, with even 
larger numerical effort by further extension to a cluster- 
RISB framework , one could introduce a two-site clus¬ 
ter in real space for each pair of Ti sites linked by an OV. 
Then possible singlet formation, contrary to the here- 
described local pair-AFM state, via intersite self-energies 
would be describable. Note that we have limited our 
work to single-OV defects. An investigation of multi-OV 
configurations up to vacancy-clustered regions on very 
large real-space lattices will surely be appreciated future 
work. In this respect, experimental information on de¬ 
fect concentrations, locations, and arrangements beyond 
currently available data is greatly needed. 

Generally, controlling the defect structure within in¬ 
terfaces of oxide heterostructures will be pivotal to even¬ 
tually engineering technological applications with des¬ 
ignated response behavior. Recent findings of room- 
temperature ferromagnetism and enhanced photocat- 
alytic performance in Ti-defected Ti02 Anatase [6l| are 
a further example for the relevance and potential of 
defect-induced oxide physics beyond the high-T^ cuprate 
paradigm. The here-documented sensitivity of the (mag¬ 


netic) electronic structure to the OV concentration could 
be useful not just for designing LAO/STO-based mag¬ 
netic order at elevated temperatures. Charge and mag¬ 
netic writing, flexible spintronic switches, sensor technol¬ 
ogy, and possible multiferroic response are only a few fur¬ 
ther optional engineering directions. The idea of creating 
atom-resolved orbital and spin polarization by controlled 
defect manipulation within a well-defined interface re¬ 
gion has so far not been translated into practicable de¬ 
vice physics. Defect control of emerging interface phases 
could be complementary to the technological potential of 
adatom-driven surface phenomena. The real-space ap¬ 
proach presented here is especially suited to simulate and 
direct such design and control of challenging correlated 
materials on a nano scale. 


ACKNOWLEDGMENTS 

We are grateful to T. Kopp, N. Pavlenko and C. 
Piefke for helpful discussions. This research was sup¬ 
ported by the Deutsche Forschungsgemeinschaft through 
FOR1346 and SFB925. Computations were performed 
at Regionales Rechenzentrum (RRZ) of the University of 
Hamburg and at the JUROPA Cluster of the Jiilich Su¬ 
percomputing Centre (JSC) under Project No. hhhOS. 


[1] A. Kalabukhov, R. Gunnarsson, J. Borjesson, E. Olsson, 
T. Claeson, and D. Winkler, Phys. Rev. B 75, 121404(R) 
(2007). 

[2] W. Siemens, G. Koster, H. Yamamoto, W. A. Harrison, 
G. Lucovsky, T. H. Geballe, D. H. A. Blank, and M. R. 
Beasley, Phys. Rev. Lett. 98, 196802 (2007). 

[3] M. Salluzzo, S. Gariglio, D. Stornaiuolo, V. Sessi, 
S. Rusponi, C. Piamonteze, G. M. DeLuca, M. Minola, 
D. Marre, A. Gadaleta, H. Brune, F. Nolting, N. B. 
Brookes, and G. Ghiringhelli, Phys. Rev. Lett. Ill, 
087204 (2013). 

[4] Z. Q. Liu, C. J. Li, W. M. Lu, X. H. Huang, Z. Huang, 
S. W. Zeng, X. P. Qiu, L. S. Huang, A. Annadi, J. S. 
Ghen, J. M. D. Goey, T. Venkatesan, and Ariando, Phys¬ 
ical Review X 3 , 021010 (2013). 

[5] A. David, Y. Tian, P. Yang, X. Gao, W. Lin, A. B. S. amd 
J.-M. Zuo, W. Prellier, and T. Wu, Scientific Reports 5, 
10255 (2015). 

[6] A. F. Santander-Syro, O. Copie, T. Kondo, F. Fortuna, 

S. Pailhe, R. Weht, X. G. Qiu, F. Bertran, A. Nicolaou, 
A. Taleb-Ibrahimi, P. L. Fevre, G. Herranz, M. Bibes, 
N. Reyren, Y. Apertet, P. Lecoeur, A. Barthelemy, and 
M. J. Rozenberg, Nature 469, 189 (2011). 

[7] W. Meevasana, P. D. C. King, R. H. He, S.-K. Mo, 
M. Hashimoto, A. Tamai, P. Songsiriritthigul, F. Baum- 
berger, and Z.-X. Shen, Nat. Mat. 10, 114 (2011). 

[8] S. M. Walker, A. de la Torre, F. Y. Bruno, A. Tamai, 

T. K. Kim, M. Hoesch, M. Shi, M. S. Bahramy, P. C. 
King, and F. Baumberger, Phys. Rev. Lett. 113 , 177601 
(2014). 

[9] W. D. Rice, P. Ambwani, M. Bombeck, J. D. Thompson, 


G. Haugstad, C. Leighton, and S. A. Crooker, Nature 
Mat. 13 , 481 (2014). 

[10] N. Reyren, S. Thiel, A. D. Caviglia, L. F. Kourkoutis, 
G. Hammerl, C. Richter, C. W. Schneider, T. Kopp, A.- 
S. Riietschi, D. Jaccard, M. Gabay, D. A. Muller, J.-M. 
Triscone, and J. Mannhart, Science 317 , 1196 (2007). 

[11] A. Brinkman, M. Huijben, M. van Zalk, J. Huijben, 
U. Zeitler, J. C. Maan, W. G. van der Wiel, G. Rijnders, 
D. H. A. Blank, and H. Hilgenkamp, Nature Mater. 6, 
493 (2007). 

[12] Y. Li, S. N. Phattalung, S. Limpijumnong, J. Kim, and 
J. Yu, Phys. Rev. B 84, 245307 (2011). 

[13] Ariando, X. Wang, G. Baskaran, Z. Q. Liu, J. Huijben, 
J. B. Yi, A. Annadi, A. R. Barman, A. Rusydi, S. Dhar, 
Y. P. Feng, J. Ding, H. Hilgenkamp, and T. Venkatesan, 
Nature Commun. 2 , 188 (2011). 

[14] J.-S. Lee, Y. W. Xie, H. K. Sato, C. Bell, Y. Hikita, H. Y. 
Hwang, and C.-C. Kao, Nature Mat. 12 , 703 (2013). 

[15] R. Pentcheva and W. E. Pickett, Phys. Rev. B 74 , 035112 
(2006). 

[16] N. Pavlenko, T. Kopp, E. Y. Tsymbal, G. A. Sawatzky, 
and J. Mannhart, Phys. Rev. B 85, 020407(R) (2012). 

[17] J. Shen, H. Lee, R. Valenti', and H. O. Jeschke, Phys. 
Rev. B 86, 195119 (2012). 

[18] C. Lin and A. A. Demkov, Phys. Rev. Lett. Ill, 217601 
(2013). 

[19] F. Lechermann, L. Boehnke, D. Grieger, and C. Piefke, 
Phys. Rev. B 90, 085125 (2014). 

[20] M. Breitschaft, V. Tinkl, N. Pavlenko, S. Paetel, 
C. Richter, J. R. Kirtley, Y. C. Liao, G. Hammerl, V. Ey- 
ert, T. Kopp, and J. Mannhart, Phys. Rev. B 81, 153414 



11 


( 2010 ). 

[21] Z. Ristic, R. DiCapua, F. Chiarella, G. M. DeLuca, 
I. Maggio-Aprile, M. Radovic, and M. Salluzzo, Phys. 
Rev. B 86, 045127 (2012). 

[22] W. Sitaputra, N. Sivadas, M. Skowronski, D. Xiao, and 

R. M. Feenstra, Phys. Rev. B 91 , 205408 (2015). 

[23] A. Joshua, J. Ruhman, S. Pecker, E. Altman, and 

S. Ilani, PNAS 110 , 9633 (2013). 

[24] K.-J. Zhou, M. Radovic, J. Schlappa, V. Strocov, R. Pri¬ 
son, J. Mesot, L. Patthey, and T. Schmitt, Phys. Rev. 
B 83, 201402(R) (2011). 

[25] J. Park, B.-G. Gho, K. D. Kim, J. Koo, H. Jang, K.-T. 
Ko, J.-H. Park, K.-B. Lee, J.-Y. Kim, D. R. Lee, C. A. 
Burns, S. S. A. Seo, and H. N. Lee, Phys. Rev. Lett. 
110 , 017401 (2013). 

[26] G. Berner, M. Sing, H. Fujiwara, A. Yasui, Y. Saitoh, 
A. Yamasaki, Y. Nishitani, A. Sekiyama, N. Pavlenko, 

T. Kopp, G. Richter, J. Mannhart, S. Suga, , and 

R. Claessen, Phys. Rev. Lett. 110 , 247601 (2013). 

[27] J. A. Sulpizio, S. Ilani, P. Irvin, and J. Levy, Annu. Rev. 
Mater. Res. 44, 117 (2014). 

[28] F. Bi, M. Huang, S. Ryu, H. Lee, C.-W. Bark, G.-B. E. 
amd P. Irvin, and J. Levy, Nature Gommun. 5, 5019 
(2014). 

[29] G. W. Bark, P. Sharma, Y. Wang, S. H. Baek, S. Lee, 

S. Ryu, G. M. Folkman, T. R. Paudel, A. Kumar, S. V. 
Kalinin, A. Sokolov, E. Y. Tsymbal, M. S. Rzchowski, 
A. Gruverman, and C. B. Eom, Nano Lett. 12, 1765 
( 2012 ). 

[30] I. R. Shein and A. L. Ivanovskii, Phys. Lett. A 371, 155 
(2007). 

[31] R. L. Johnson-Wilke, D. Marincel, S. Zhu, M. P. Waru- 
sawithana, A. Halt, J. Sayre, K. T. Delaney, R. Engel- 
Herbert, C. M. Schlepiitz, J.-W. Kim, V. Gopalan, N. A. 
Spaldin, D. G. Schlom, P. J. Ryan, and S. Trolier- 
McKinstry, Phys. Rev. B 88, 174101 (2013). 

[32] J. G. Li, J. 1. Beltran, and M. C. Munoz, Phys. Rev. B 
87 , 075411 (2013). 

[33] A. Janotti, J. B. Varley, M. Ghoi, and G. G. Van de 
Walle, Phys. Rev. B 90, 085202 (2014). 

[34] A. Lopez-Bezanilla, P. Ganesh, and P. B. Littlewood, 
Phys. Rev. B 92, 115112 (2015). 

[35] D. Grieger and F. Lechermann, Phys. Rev. B 90 , 115115 
(2014). 

[36] K. Michaeli, A. C. Potter, and P. A. Lee, Phys. Rev. 
Lett. 108 , 117003 (2012). 

[37] G. Chen and L. Balents, Phys. Rev. Lett. 110 , 206401 
(2013). 


[38] S. Banerjee, O. Erten, and M. Randeria, Nature Phys. 
9, 626 (2013). 

[39] J. Ruhman, A. Joshua, S. Ilani, and E. Altman, Phys. 
Rev. B 90 , 125123 (2013). 

[40] L. Yu and A. Zunger, Nat. Common. 5, 5118 (2014). 

[41] C. Lin and A. A. Demkov, Phys. Rev. Lett. 113 , 157602 
(2014). 

[42] N. Pavlenko, T. Kopp, and J. Mannhart, Phys. Rev. B 
88, 201104(R) (2013). 

[43] C. Zener, Phys. Rev. 82, 403 (1951). 

[44] P. W. Anderson and H. Hasegawa, Phys. Rev. 100, 675 
(1955). 

[45] H. Y. Hwang, Y. Iwasa, M. Kawasaki, B. Keimer, N. Na- 
gaosa, and Y. Tokura, Nature Materials 11 , 103 (2012). 

[46] A. N. Rubtsov, V. V. Savkin, and A. 1. Lichtenstein, 
Phys. Rev. B 72 , 035122 (2005). 

[47] P. Werner, A. Comanac, L. de’ Medici, M. Troyer, and 
A. J. Millis, Phys. Rev. Lett. 97, 076405 (2006). 

[48] O. Parcollet, M. Ferrero, T. Ayral, H. Hafermann, 
1. Krivenko, L. Messio, and P. Seth, Computer Physics 
Communications 196 , 398 (2015). 

[49] L. Boehnke, H. Hafermann, M. Ferrero, F. Lechermann, 
and O. Parcollet, Phys. Rev. B 84, 075145 (2011). 

[50] F. Lechermann, A. Georges, G. Kotliar, and O. Parcol¬ 
let, Phys. Rev. B 76 , 155102 (2007). 

[51] B. Amadou, F. Lechermann, A. Georges, F. Joliet, T. O. 
Wehling, and A. 1 . Lichtenstein, Phys. Rev. B 77 , 205112 
(2008). 

[52] T. Li, P. Wolfle, and P. J. Hirschfeld, Phys. Rev. B 40, 
6817 (1989). 

[53] J. Biinemann, W. Weber, and F. Gebhard, Phys. Rev. 
B 57 , 6896 (1998). 

[54] F. Lechermann, Phys. Rev. Lett. 102 , 046403 (2009). 

[55] L. Huang, L. Du, and X. Dai, Phys. Rev. B 86, 035150 

( 2012 ). 

[56] E. C. Andrade, E. Miranda, and V. Dobrosavljevic, 
Phys. Rev. Lett. 102 , 206403 (2009). 

[57] G. Kotliar and A. E. Ruckenstein, Phys. Rev. Lett. 57, 
1362 (1986). 

[58] F. Lechermann, S. Biermann, and A. Georges, Prog. 
Theor. Phys. Suppl. 160 , 233 (2005). 

[59] M. R. Fitzsimmons, N. W. Hengartner, S. Singh, M. Zh- 
ernenkov, F. Y. Bruno, J. Santamaria, A. Brinkman, 
M. Huijben, H. J. A. Molegraaf, J. de la Venta, , and 
1. K. Schuller, Phys. Rev. Lett. 107 , 217201 (2011). 

[60] A. Ron, E. Maniv, D. Graf, J.-H. Park, and Y. Dagan, 
Phys. Rev. Lett. 113 , 216801 (2014). 

[61] S. Wang, L. Pan, J.-J. Song, W. Mi, J.-J. Zou, L. Wang, 
and X. Zhang, J. Am. Chem. Soc. (2015). 


